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ABSTRACT 


The  acoustic  radiation  characteristics  of  a  fluid  filled  concentric  shell 
structure  are  investigated  by  making  use  of  the  simple  source  surface  integral 
method.  The  fluid  field  is  described  in  terms  of  surface  expansion  functions, 
while  the  in  vacuo  structural  modes,  obtained  by  using  the  B0S0R4  structural 
program,  are  used  as  a  basis  for  the  shell  motions.  The  structural  code  can  be 
applied  to  branched  shells  of  revolution,  thereby  obviating  the  need  to 
introduce  structural  coupling  of  the  two  shells  in  a  manner  which  is  external 
to  the  main  code . 

For  illustrative  purposes  three  types  of  internal  harmonic  excitations 
have  been  considered:  (1)  an  axi symmetric  unit  pressure  acting  on  the  central 
region,  (2)  a  unit  line  load  acting  on  the  central  region,  and  (3) 
concentrated  axial  loads  acting  at  the  poles  of  the  spherical  closures.  The 
far  field  fluid  pressures  have  been  obtained  for  a  wide  range  of  frequencies 

for  both  the  single  shell  and  for  the  double  shell  with  and  without  inner 
fluid. 

Because  of  numerical  difficulties  encountered  when  applying  the  integral 
methods  to  internal  fluid  problems,  some  simple  approximations  to  the  internal 
fluid  problem  were  considered.  The  results  obtained  compare  quite  favorably 
with  those  of  the  exact  solution  over  a  wide  range  of  frequencies. 
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LIST  OF  SYMBOLS 

a  (b) 

radius  of  inner  (outer) shell 

Aj  (A^) 

surface  area  of  inner  (outer)  shell 

c 

acoustic  velocity  in  fluid 

d 

distance  between  field  and  source  points 

pressure  influence  coefficient  associated  with  k  th 
field  point  and  j  th  source  point 

velocity  influence- coefficient  associated  with  k  th 
field  point  and  i  th  source  point 

G 

free  space  Green's  function 

G 

n 

n  th  circumferential  component  of  Green’s  function 

h 

distance  between  inner  and  outer  shells 

k 

wave  number 

K 

generalized  stiffness  of  the  structure 

1 

length  of  cylindrical  region 

m,  n 

axial  and  circumferential  mode  numbers,  respectively 

M 

generalized  mass  of  the  structure 

n 

P 

normal  to  surface  at  location  P 

P 

fluid  pressure 

q 

q 

generalized  coordinate 

generalized  coordinate  associated  with  surface  expansion 
function 

Q 

generalized  force 

Q 

generalized  force  associated  with  surface  expansion  function 

r,  0,  X 

radial,  circumferential  and  axial  coordinates,  respectively 

R 

mean  radius  (R  =  ^  ^ 

s 

meridional  arc  length  of  surface  of  revolution 
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w 


m 


V 

p 

'p 

03 

0) 

m 


source  strength  distribution  function 

axial »  circumferential  and  normal  components  of  ni  th 
structural  mode 

normal  displacement  component  of  shell 

impedance  function  associated  with  j  th  surface  expansion 
function 

radial  coordinate  measured  from  mean  radius  R 

Laplacian  operator 

fluid  density 

velocity  potential 

k  th  surface  expansion  function 

frequency  of  excitation 

in-vacuo  frequency  of  m  th  structural  mode 


I.  INTRODUCTION 


A  vast  amount  of  literature  exists  concerned  with  the  prob¬ 
lem  of  predicting  sound  radiation  from  submerged  structures  exe¬ 
cuting  induced  harmonic  vibrations. 

Except  for  some  special  geometries,  for  which  analytical 
solutions  can  be  achieved  by  the  standard  technique  of  separation 
of  variables,  the  great  majority  of  investigations  have  been 

out  by  using  integral  equation  methods.  Although  some 
solutions  have  been  obtained  by  applying  the  finite  element 
methods  these  methods  do  not  compare  favorably  with  the  simplici¬ 
ty  and  economy  resulting  from  the  use  of  the  integral  methods, 
especially  so,  because  the  problems  of  interest  are  generally 
those  for  which  the  fluid  field  is  of  infinite  domain.  The  form 
of  the  integral  equation  method  most  often  used  to  solve  the 
radiation  problem  is  referred  to  as  the  "simple  source  method". 
Here  one  assumes  that  the  acoustic  field  potential  is  expressed 
in  terms  of  a  layer  of  monopole  sources  distributed  over  the 
radiating  surface. 

Upon  invoking  the  compatibility  relationship  for  the  normal 
velocities  at  the  interface,  one  obtains  an  integral  equation  for 
the  unknown  source  density  function.  Once  the  density  function 
is  determined,  one  can  compute  the  pressures  and  velocities 
throughout  the  fluid. 

It  should  be  noted  that  the  use  of  the  integral  equation 
methods  does  present  difficulties  at  and  in  the  neighborhood  of 
the  characteristic  frequencies  associated  with  the  interior 
Dirichlet  problem.  However,  it  has  been  pointed  out  by  Piacyzyk 
and  Klosner  [1]  that  the  internal  resonance  problem  which  is 


Indeed  present  In  the  transducer  type  problea;  l.e.,  when  the 
velocity  of  the  radiating  surface  is  specified,  does  not  occur  In 
the  structural  radiation  problea  If  the  coupled  equations  are 
solved  in  the  proper  order. 

In  almost  all  of  the  Investigations  of  structural  radiation 
the  fluid  is  In  contact  with  only  one  side  of  the  radiating 
surface.  There  is  ,  however,  a  paucity  of  literature  concerned 
with  radiation  froa  structures  having  Internal  and  external 
wetted  surfaces  which  Interact  with  one  another.  One  such 
investigation,  carried  out  by  Cambridge  Acoustical  Associates 
considered  the  radiation  from  a  structure  consisting  of  two  flat 
plates  separated  by  longitudinal  structural  membranes,  and  filled 
with  a  fluid.  One  side  of  the  composite  structure  was  in  contact 
with  a  semi-infinite  fluid,  while  the  other  side  was  in  vacuo. 

A  preliminary  investigation  of  a  more  realistic  situation  was 
conducted  aat  the  Polytechnic  Institute  of  New  York.  They 
considered  the  acoustic  radiation  from  a  structure  consisting  of 
concentric  spherically  capped  cylindrical  shells  separated  by, 
and  immersed  in,  an  acoustic  fluid.  The  Internal  shell  was 
stiffened,  and  the  two  shells  were  connected  both  structurally  and 
acoustically. 

In  the  present  investigation  the  problem  of  radiation  from  a 
concentric  shell  structure  is  reconsidered  by  using  a  more 
refined  and  efficient  approach. 

First,  surface  expansion  functions  are  used  to  describe  the 
fluid  fields.  The  expediency  of  using  expansion  functions  have 
been  amply  demonstrated  for  both  shock  [2]  and  radiation  [3] 
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problems.  Their  use  results  in  reduced  core  memory  requirements, 
thereby  allowing  for  a  more  detailed  description  of  the  fluid 
fields,  which  is  so  necessary  in  order  to  accurately  describe,  in 
particular,  the  azisymmetric  response  of  the  inner  fluid. 

Second,  the  in  vacuo  structural  modes,  which  are  used  to 
describe  the  shell  motions,  are  obtained  by  using  the  BOSOR  4 
structural  program.  This  highly  efficient  program  uses  a  circum¬ 
ferential  modal  expansion,  and  can  be  applied  to  branched  shells 
of  revolution,  thereby  eliminating  the  need  to  introduce  struc- 
tural  coupling  of  the  two  shells  in  a  manner  which  is  external  to 
the  main  code. 


II.  inmEX 


I  Consider  the  acoustic  radiation  characteristics  of  a  shell 
submerged  in  an  acoustic  medium  and  subjected  to  a  time-harmonic 
excitation. 

The  velocities  at  the  fluid-structural  interface  are 
expressed  by  sets  of  orthogonal  functions.  For  the  fluid, 
harmonic  surface  expansion  functions  are  used,  while  for  the 
structure  the  in  vacuo  structural  modes  serve  as  an  expansion 
basis. 

In  the  following,  we  will  describe  and  formulate  the 
procedure. 


II. 1  Fluid.  Field 

The  Helmholtz  wave  equation  for  an  inviscid  fluid 
undergoing  simple  harmonic  motion  is 

(V^  +  k^)4)  -  0 

where  (|)  =  fluid  velocity  potential;  the  wave  number 

k  =  u/c;  c=acoustic  velocity. 

The  particle  velocity  is  expressed  as  the  gradient  of  <|> 

V  *  -V4) 

and  the  pressure 


(2-1) 


(2-2) 

(2-3) 


The  solution  of  Eq.  (2-1)  corresponding  to  a  unit 
source,  i.e,  the  free  space  Green's  function  representing  a 
divergent  spherical  wave,  can  be  written  as 

P  -Ikd 

b  *  -  e _ 

4iTd 

where  d=distance  from  source  point  Q  to  field  point  P.  (Fig.  1) 

For  convenience,  the  time  factor  is  supressed  here  and 

is  omitted  in  what  follows,  unless  clarity  demands  its  inclusion. 


(2-4) 
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By  invoking  the  Helmholtz  integral 


e<{>(P) 


ii'll 


1-^  d  al 
a  dn 


(2-5) 


where  e  »  0,  for  P  not  V,  e  =  J5  for  P  on  boundary  a,  and  £  ■  1,  P  in  V, 
we  can  express  the  value  of  tp  at  a  point  P  in  the  region  V  in 

terms  of  (}>  and  its  normal  derivative  34)/3n  at  the  fluid  bound¬ 
ing  surface  a  .  The  first  integral  of  Eq,  (2-5)  represents  the 
contribution  to  the  potential  due  to  a  distribution  of  simple 
sources  on  the  bounding  surface,  while  the  second  term  represents  that 
due  to  distributed  doublets  wi*th  axes  normal  to  the  surface. 

We  can,  under  certain  conditions,  express  the 
value  of  the  potential  in  the  region  in  terms  of  simple  sources 
only,  or  doublets  only,  distributed  over  the  fluid  boundary. 


Expressed  in  terms  of  simple  sources,  the  potential 


^(P)  -  - 


S  -ikd  . 
:r  e  da 

a 


(2-6) 


where  S  denotes  the  source  strength  distribution  function.  Note 
that  when  point  P  =  p  lies  on  the  surface  the  integrand  is 
singular,  and  the  normal  derivative  becomes 

-ikd 


3n 


(p)  *  j  S(p)  -  PV 


J 


^  3n  ''  d 


-)da  1 


(2-7) 


where  PV  =  principal  value. 


When  considering  surfaces  of  revolution,  (Fig.  2), 
it  is  expedient  to  use  a  circumferential  harmonic  expansion,  so  that 


*This  expression  is  invalid  when  k  corresponds  to  the 
internal  Dirichlet  or  Neumann  eigenvalue  problems  (See  Ref.  4  , 

I 

pp.  499,  also  Ref  1  ) 
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S  -  ^  S  COi  n3 
n  n 


y 

0  -  „  '*’n  C04  ne 


where 


s 


n 


G^(P,Q)S^(Q)ds 
n  n 


and  the  pressure  influence  coefficient 


2u  -ikd 
e 


a06  nJp  dJp 


o  d 

s  is  the  arc  length  measured  along  the  generator;  r(Q)  =  distance 
from  axis  to  generator  at  source  point  Q.  Eq.  (2-7)  is  now 
written  as 


3()), 


n(p)  _  1 


s 


Y  s  (p)  +  PV  [  —  [G^(p,Q)]S^(Q)ds] 

•^o  p 


and  the  velocity  influence  coefficient. 


[G„(P.Q)]-  -  r(Q) 


9n  ^  n'  , 

p  477 


277 


-ikd 


—  )  C04  n0  d9 


^  o  p 

Of  course^  if  the  field  point  P  does  not  lie  on  the  surface,  then 


94)n  (P) 
9n 


j  3r^  n 


[G„(P,Q)]S„(Q)ds 


n 


Upon  subdividing  the  surface  into  N  bands  along 
the  generator,  and  assuming  that  the  source  strength  is  constant 

I 

along  the  generator  within  each  band,  we  obtain  the  discretized 
forms  of  Eqs.  (2-10),  (2-12)  and  (2-14) 

I  D,  ,  S  (j) 

K  j  =  l 

where  Dj^j  and  are  the  pressure  and  velocity  influence 


(2-8) 

(2-9) 

(2-10) 

(2-11) 

(2-12) 

(2-13) 

(2-14) 

(2-15a) 


(2-15b) 
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coefficients  associated  with  reference  to  a  field  point  P  in  the 
kth  band  and  reference  to  a  source  point  Q  in  the  jth  band. 

II. 2  SURFACE  EXPANSION  FUNCTIONS  -  IMPEDANCE  FUNCTIONS 


Since  we  are  concerned  with  surfaces  of 
revolution,  a  circumferential  modal  expansion  is  used.  The 
solution  for  each  circumferential  wave  number  is  obtained  and  the 
final  result  obtained  by  superposition. 

Thus  for  each  circumferential  mode,  an  orthogonal 
set  of  fluid  surface  expansion  functions  is  chosen  for  use  as  an 
expansion  basis  for  the  normal  component  of  the  surface 
displacement  w, 

w  = 


I 


=1 


iut 


(2-16) 


It  is  to  be  understood  that  all  of  the  above 
have  the  same  harmonic  circumferential  variation, and  that  they 
have  been  normalized  with  respect  to  the  surface  area  A,  so  that 


\  da  -  A 

The  impedance  function  zj  is  defined  as  the 
surface  pressure  field  p  j ,  due  to  the  surface  velocity 
is  determined  by  equating  the  normal  components  of  the  fluid  'velocity 
at  the  surface  to  that  associated  with  the  function  ip. 

By  appealing  to  Eqs.  (2-15),  we  can  express  the  above  and  the 
corresponding  impedance  function  in  matrix  form. 


(2-17) 


(2-18) 


-  [OpQ.n^HS^} 


(2-19a) 
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up  -  i-^tDpQUs^} 

The  simultaneous  solution  of  Eqs.  (2-14)  then 

yields  the  source  and  impedance  vectors  and  which  are  due 
to  the  surface  velocity  distribution- . 


In  turn,  we  represent  in  terms  of  the  surface 
expansion  functions,  so  that 


(2-19b) 


=  f  Y 

k=l 


where 


=  J 


a 


Z^ip.  da  =  n(l+6  ) 

k  on 


6  = 

on  ^0,n?^0 


fs 


r(s)Z^iijj^ds 


(2-20) 

(2-21) 

(2-22) 


Observe  that  is  the  generalized  force  associated  with  the 

-displacement  and  the  pressure  field  arising  from  \p^  . 

II. 3  STRUCTURAL-FLUID  INTERACTION 
The  normal  component  of  the  displacement  of  the  shell  can 
be  expressed  as 

w  =  y  q  w  (2-23) 

^  mm 
m=l 

where  Wjj^  =  the  w-component  of  the  mth  in  vacuo  shell  mode,  and 
q^jj  =  mth  generalized  coordinate.  Again  it  should  be  emphasized 
that  all  of  the  above  modes  are  associated  with  one  particular 
circumferential  wave  number. 

For  harmonic  time-dependent  steady-state  motion,  the  shell 
equation  of  motion  becomes 

^•m  (  Hn  ~  q_  “  Q  ^  +  0 
mm  ^  ^m 


(2-24) 
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where'^jjjs  in  vacuo  frequency  of  the  mode,  the  mth 
component  of  the  generalized  mass  and  force  is,  respectively. 


M 

m 


ph(U  +V  +W  )  da 


m 


m 


m 


(2-25) 


Q 


m 


-  P(J>  da 
JJ  ^ 
a 


(2-26) 


and  QEjjj  =  generalized  force  associated  with  the  applied  loads. 
UmrVm/Wm  are  the  shell  displacements  along  the  generator  and  in 
the  circumferential  and  normal  directions,  respectively. 

At  the  shell-fluid  interface,  the  normal  shell  velocity 
must  be  equal  to  the  normal  fluid  particle  velocity.  Hence 


M 

iuw  =  ico  ^ 

m=l 


m  m* 


K 


I 


k-1 


k^k 


(2-27) 


Making  use  of  the  orthogonality  relationship,  we  have 


where 


M 


q,  =  ^  X  q 
K  ^  1  km  m 
111=  i 


=  T  fj  "  r, 

km  A  jJ  m  k 
a 


W  d  a 


(2-28) 


(2-29) 


Expressing  the  surface  pressure  in  the  form  of  series 
expansions  in  both  the  surface  expansion  functions  and  the  in 
vacuo  structural  modes,  we  obtain  the  following 


if  Qn.  K  Q, 

p-i  x* 

m=l  ^m  k=l 


(2-30) 


where  the  associated  generalized  forces 


‘^m  =- 


p  W  da 
m 


(2-31) 
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Jj 


P  dcr 


(2-32) 


are  related  as  follows: 

k*l  ^ 


’‘I  \nk  Qi 


(2-33) 


k“l 


Recall  that  pj  denotes  the  surface  pressure  due  to  the 
surface  velocity  .  Thus,  upon  observing  Eq,  (2-16),  we 
note  that  the  presssure  due  to  the  actual  surface  velocity  can  be 
written  as  K 


P  -  I  q.  ZJ 
j=*l  J 


(2-34) 


Upon  substituting  Eq.  (2-34)  into  Eq.  (2-32),  we  obtain 


K 

-I 


q. 


d  o 


which,  upon  substituting  (2-21)  and  (2-28) ,  becomes 

M 

^  \ 

m-1 


^  K  ^  K  M 


Combining  Eqs.  (2-33)  and  (2-36)  then  yields 

K  K  M 

®  k»l  j=l  m«l  J  ‘"J 


(2-35) 


(2-36) 


(2-37) 


The  substitution  of  Eq.  (2-37)  into  the  equations  of  motion  (2.- 
24),  leads  to  the  following  set  of  coupled  simultaneous 
equations,  expressed  in  matrix  notation, 

[K  -  ti/M  'ff  L]{q}  =  {Q^} 

[K]  -  M  ] 
m  m 

[  dfM]  =  [u^M  ] 
m 


(2-38) 


where 


(2-39) 
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The  solution  of  the  above  equation  yields  the  generalized  coordinate 
vector  q  ,  By  substituting  q  into  Eq.  (2-28)  and  then  into  Eq.  (2-34) ^  we  obtain 
an  expansion  for  the  surface  pressures  at  the  i  nodes  (or  bands)  (i  =  l-^I) 


where  Z^j  =  2ij  =  the  value  of  the  impedance  function  at  the  ith 
node  associated  with  the  surface  displacement. 

The  far-field  pressure  is  then  obtained  by  substituting  the 
source  strength  vector  SJ  into  the  discretized  form  of  Eq.  (2-10) 
and  thereupon  substituting  into  Eq  (2-3).  Thus  the  far-field 
pressure  at  point  P, 

P(P)  -  ia)p{G  (P)}^[s  ][A,  ]{q  } 

^  ij  jm  m 


where 


^i(P)  =  pressure  influence  coefficient  vector,  [G^(p)  -  G(P,i)] 
^ij  =  element  of  source  strength  vector 

4 

As^  =  product  of  ith  band  width  and  value  of  at  ith  node] 
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III. APPLICATION  TO  THE  DOUBLE  HULL  SHELL 

We  now  consider  the  concentric  capped  cylindrical 
shell  (Fig  3),  and  define  a  set  of  surface  expansion  functions 
for  the  inner  and  outer  shells  and  ,  respectively. 

K>  1C 

Following  along  the  lines  described  in  Section  ll,  we  have  that 


I 


w 


M 


I 

m=l 


q 

m  ni 


TT  « 

w“=  I 

m 


II 


^  m 


K2 


I 

k=l 


(3-1) 


(3-2) 


where  the  subscript  I  and  II  refers  to  the  inner  and  outer 
shell,  respectively,  and  all  of  the  above  functions  have  a 

common  circumferential  wave  number  n.  Wm  is  the  normal  component 
of  the  mth  in  vacuo  mode  of  the  branched  shell  structure,  and 

m 

and  refer  to  that  part  of  Wm  associated  with  the  surfaces 

of  the  inner  and  outer  shells,  respectively.  Kl  and  K2  are  the 
number  of  expansion  functions  used  for  the  inner  and  outer 
surfaces,  respectively. 

As  in  Section  II,  we  have  that 

1-1  1 
K  ^  ,  km  m 
m=l 


(3-3) 


(3-4) 


where 


km 


1 

A, 


da 


m 


(3-5) 


2^ 


m  k 


(3-6) 
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We  define  the  impedance  functions  as  the  pressure  at  the 

kth  surface  due  to  the  surface  expansion  at  the  1-th  shell 

surface.  It  is  evaluated  by  following  a  procedure  similar  to 

that  for  a  single  shell*.  For  the  inner  fluid  the  discretized 

form  of  the  compatibility  equation  (2-12)  is  applied  to  all 

wetted  nodes  of  the  inner  fluid;  for  example,  when  the  fluid 

velocity  is  it  is  implied  that  e  o  at  all  nodal 

j  dn 

points  lying  on  the  outer  shell.  When  considering  the 
outer  fluid,  the  impedance  function  is  determined  by  using  the 
procedure  described  for  the  single  shell.  Thus  for  the  inner 
fluid  the  impedance  functions  associated  with  the  displacement  4^^ 
along  the  inner  shell  are  expressed  as 


11 


K1 

I 

k=l 


11 

Zki 


(3-7) 


(3-8) 


*As  a  consequence  of  the  fact  that  the  normal  direction  is  always 
assumed  to  be  positive  outward,  the  term  1/2  S(p)  which  appears 
in  Eq.  (2-7)  is  replaced  by  -1/2  S(p)  when  considering  the  inner 
fluid  in  contact  with  the  outer  shell. 
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while  those  associated  with  the  displacement  distribution 


II 


of  the  outer  shell  are 


22 

K2  yf.  , II 

22  ^  A 

k-1  ^2 


(3-9) 


where 


K1  Y 


12 

kj  .  I 


'12 


k»l  ^1 


y 


11 


kj  Jj 


4i  fk  '“s 


(3-10) 


(3-11) 


1 


(3-12) 


and 


22 

\j 


JJ 


^2 


(3-13) 


Y 


.12 

kj 


i  ^^2  < 


(3-14) 


The  subscripts  and  superscripts  1  and  2  respectively  denote  the 
surfaces  of  the  inner  and  outer  shells  which  are  in  contact  with 
the  inner  fluid. 

Finally,  for  the  outer  fluid 


(3-15) 


(3-16) 


One  can  readily  show  that  the  components  of  the  generalized  force 
associated  with  the  interaction  fluid  pressures  are 
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K1  K1  M  ,, 

Q  - -I  I  I  ^ 

U  j-l  m-1  ^  ” 

K2  K1  M  tt  91  T  ^ 

“  k-l  i-l  il  “ 


.22 


K2  K2  M 

^  ^  ^  ,  'kj  “jm  m 

k=l  J»1  m-1  nk 


K1  K2  M  _  , 

'Ll  j-i  Li  "j”  “ 


.12  ^II  q 


(3-17) 

(3-18) 

(3-19) 

(3-20) 


and 


K2  k: 

-I  I 


k-l  j=l  Ll 


ill  -.33  ,II 


(3-21) 


where  the  negative  signs  appear  on  the  right  hand  sides  of  the 

11  12  33 

expressions  for  ,q^  ,q^  because  positive  pressures  are  in  the 
opposite  direction  to  positive  radial  displacements  for  these  cases. 
The  shell  equation  for  time-harmonic  motion  is 


M  (i/ 
m  m 


- 


q  .  Q*  +,,11+  q21  q.  q22  q.  ql2  *  q33 

m  m  m  m  m  m  m 


[  K-  0)  M  -  L]  {q}  -  {Q^} 


(3-22) 

(3-23) 


A 


where  K  and  U)  w  have  been  previously  defined  [Eq  2-39  ]  and 
here 

utilizing  the  q-matrix  found  from  Eq.  (3-23)  ,  we  can 
determine  the  pressure  at  the  surface  nodes  from  the  following 
expressions 


(3-24) 


(3-25) 
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.nd  {p”)  -  [ZJ33  -  zJj]  »-26) 

The  far-field  pressures  are  obtained  by  simply  applying  the  form 
given  by  Eq.  (2-41.)  to  the  outer  fluid. 
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IV. APPROXIMATE  SOLUTIONS 

Serious  numerical  difficulties  are  encountered  when  applying 
the  integral  methods  to  solve  radiation  problems  for  an  enclosed 
fluid.  As  a  consequence,  the  number  of  nodes  (or  bands)  required 
to  obtain  accurate  results  must  be  greatly  increased  over  that 
which  is  required  for  the  external  fluid.  In  addition  much 
greater  detail  has  to  be  introduced  when  evaluating  the  fluid 
influence  coefficient  integrals.  Because  of  this,  a  preliminary 
investigation  of  some  simple  approximations  was  undertaken. 


Consider  a  fluid  region  within  two  concentric 
cylindrical  surfaces  and  end  closure  rings  as  shown  in  Fig.  4. 
Let  the  normal  displacement  of  the  inner  surface  be  of  the  form 
cos  cos  ,  and  let  the  outer  surface  and  the  end 

closure  rings  be  rigid.  For  this  problem  first  the  exact 
solution  is  obtained  and  then  approximations  of  the  fluid  field 
are  developed  for  the  situation  when  the  gap  between  the  two 
cylinders  is  small. 

IV. 1  Exact  SoIutlQil 

The  fluid  field  potential  is  represented  by  the 


wave  equation 


2  2 
O  4. 

2  ^  2 
99^ 


(4-1) 


where  r,  9,  x,  are  cylindrical  coordinates.  For  the  problem 
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described  above,  the  boundary  conditions  are. 


8((>^  \ 

»  -i(jJ  C04  (^^)C04  n9 

(4-2) 

If(r-b)  -  0 

(4-3) 

1^  (x-o,-£)  -  0 

(4-4) 

The  general  solution  of  Eq.  (4-1),  satisfying  the 
boundary  condition  in  Eq.  (4-4)  is 


(4-5) 


where 


k 

m 


2  2  2  2  2 
m  IT  j.  0)  ^  m  TT 

-2“  ~r-~i 

9^  9^ 


< 


(4-6a) 


^nr  Sn  are  Bessel  functions  of  the  1st  and  2nd  kind  of  order 
n. 


or 


for 


2  2 
m  7T 


(4-6b) 


Sn  are  the  modified  Bessel  functions. 

The  constants  Amn  and  Bmn  determined  from  Eqs.  (4-2)  and 
(4-3) ,  and  thus  Eq.  (4-5)  can  be  written  as 


=  I  I 


iO) 


m=0  n=0  m 


(k  b)g  (k  r)  g  (k  b)f  (k  r) ,  .  //  -7^ 

[  n  m  m  -  n  m  n  m  ]  ^  ,m7TXx  i6L5t  (4-7) 

- - - i  cos  n9  cos  (-y“)e 

nm  nm  nm  nm* 
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where 


9f„ 
3(  ) 


The  velocity  and  pressure  fields  are  then  obtained  by 
substituting  into  Eqs.  (2-2)  and  (2-3) 

IV. 2  FIRST-ORDER  APPROXIMATION  -  Series  Expansion 

For  the  region  shown  in  Fig.  4,  we  express  the  fluid 
potential  in  terms  of  a  three-term  series  expansion  through  the 
thickness , 

<P  *  f^(x,e,t)+  zf j^(x,e,t)+z^f2(x,e,t) 


where 


the  mean  radius 


and  the  gap 


z  -  r-R(-h  <  z  ^  h  ) 
2  -  -  2 


R 


a+b 


h  »  b-a 


(4-8) 

(4-9) 

(4-10) 

(4-11) 


The  fluid  equations  of  motion  can  be  expressed  in 
terms  of  Hamilton's  variational  principle. 


5  [  I  “  Ldt  + 

t. 


Wdt  ]  =0 


(4-12) 


'1 


where  the  Lagrangian 

L  -  j||tPl(|^)'  +  (He)'  +  (f)'l  -  ^ 

and  the  work  done  by  the  pressure  at  the  boundaries 


(4-13) 


w.  ±  I 


n  ^ 

IF 


(4-14) 


The  (+)  sign  refers  to  the  surfaces  r=a  and  x=o,  and  the  (-)  sign 
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refers  to  the  surfaces  r=b  and  x  -1.,  and  s  normal  surface 
velocity  (positive  in  positive  r  and  x  directions). 

The  variational  equation  (4-12)  can  thus  be  written  as 


|i)rdedrdx,dt  i  P 


)) 


It  ^ 

.  (4-15) 


which  upon  integrating  by  parts  and  enforcing  the  condition  that 
the  variation  vanishes  at  the  end  points,  becomes 

-Pf'"  ^  +  ^4)  5*rdedrdx]dt 


2  „  2  2-2 
r  99  dx  c  9t 


+  P 


[//  (f  ^ 


z=h/2 

z=-h/2 


rdGdx  +  1  I  (-1^  +  V  )6(}) 
j )  ax  n 


=£ 


X' 


x=o 


rd6dr]dt  =  0 


(4-16) 


The  evanescence  of  the  above  surface  integral  simply  leads  to 
the  interface  boundary  conditions 


and 


94) 

9r 


94) 

9x 


V  at  z=±  h/2 
n 


-V  at  x=o,-£.  . 
n 


(4-17) 


If  we  introduce  Eq.  (4-8)  into  Eq.  (A-l?)'  ,  then 
-v^(b)  (i).  II  (h/2)  -  (^  +  hfj 

■  H  "If  ■  *'^2 


and  thus 


[v^(b)  +  v^(s)l 

[v  (b)  -  V  (a)] 
n  n 


(4-18) 


(4-19) 
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Since  f2  are  expressed  in  terms  of  the  surface 

displacements  t 

<S4>  “ 


(4-20) 


It  should  be  noted  that  this  procedure  can  easily  be  generalized 
to  include  a  higher  -  order  expansion  function. 

Upon  utilizing  the  approximation 

'  1_  ~  i  ^ 

R+z  ■  R  " 

and  integrating  through  the  thickness,  Eq.  (4-16)  takes  the  following  form 

t„ 


(4-21) 


[ 


'1 


{2Rhf2  ll  99  _  2i_  2^1  604.’  —  ^2  69  +  h^  f,  + 

Z  i  ^  o.OO  -^2r^1,00+  Z,H9  +  o.xx  — 


’l,xx 


IT  *2-’“  "t  '“'o.tt  +  4  'l.tt  ^  ‘  ° 

c 


(4-22) 


Setting  the  above  integrand  equal  to  zero,  we  obtain  the 
differential  equation  for  fo' 

o,xx  ■  “2  ^o,tt  “  "  ^^2  ■  R  ^1  "l2Rf"  R  ^2,66  ■’’•^l.xx'^ 

R  R 


+  Rf„  -  1  f,  -R  f_ 

2, XX  “2  l.tt  -y  2,tt] 

c  c 


(4-23) 


Since  h/R  <<1,  the  terms  appearing  in  the  bracket  are 
neglected  and  the  final  equation  is 


f  +  1  f  00  1  f  =  -2f,  -  i-  f, 

0»XX  “;r  O»00  -  — ;r  0,tt  2  R  1 


Letting 


(4-24) 


iwt 


f^  =  f^  (x)  cos  n6  e 


(4-25) 
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we  have 


o«xx 


2  2 


c2 


(4-26) 


We  can  reformulate  the  problem  in  terms  of  the  pressure,  so  that 


P  “  (p^  +  z  pj^  +  P2)  cos  n9 

and  the  equivalent  of  Eqs.  (4-26)  and  (4-19)  are 

2  2 

‘^o.xx  ^2  ^^2*^0  z  ^  t-i 


(4-27) 


(4-28) 


and 


p'  »  -ipu  [^n(b)  +  ^n(a)  ] 


p  -  -ipw  [  n(b)  -  n(a)] 

^  2h 


(4-29) 


For  the  problem  described  in  the  beginning  of  this 


section  (Fig.  4)  ,  ^  /-n  .  tucos  m-rrx  cos  n9 

n''  -JT 


V  (b)  =  0 
n 


(4-30) 


and  since  the  end  closure  rings  are  assumed  to  be  rigid 


or  alternatively 


-|i*0atx*o,  t 
8x 


-  0  at  X  -  o,  ^ 
dx 


(4-31) 

(4-32) 


By  substituting  Eqs.  (4-30)  into  (4-29),  we  note  that  and  P2 
identically  satisfy  this  last  condition,  and  it  is  thus  necessary 
that 


3x 


0  at  X  «  o,  £ 


(4-33) 


Applying  these  end  conditions  to  the  solution  of  Eq.  (4-28), 
results  in  the  following 


Po 


{ 


2apa) 


»,2  2 
b  -  a 


/  2  2 
/  rJiL  _  IL.  _ 


(T“)  J  }  cos 


mTTx 


(4-34) 


-23- 


The  pressures  can  thereupon  be  expressed  as 


2a 


p  ”  pcu  {  2  2  2  2 

[b  -  a  ]  .0)  n 


r  uj  11  ,ni7r  X  2 1 

-  ::2  ■  ' 
c  R 


+  1  Z  _  z _ 

^2  2(b-a) 


}  cos  mTTx  cos  ,, 

(4-35) 


IV.  3  ONE-DIMENSIONAL  APPROXIMATION  -  "PIPE  FLOU" 


The  linearized  form  of  the  continuity  and  momentum  equations 
for  inviscid,  compressible  flow  are: 


+  p  V  •  V.  0 

(4-36) 

-  i  Vp  .  ^ 

p  9t 

(4-37) 

Since  the  acoustic  velocity 

c  =(i2_)l/2 
'■9p  ^ 

(4-38) 

we  can  rewrite  the  continuity  equation  as 

\  l^+P  V.v  =  o. 


where  in  cylindrical  coordinates 


V 


V 


J. 

r 


■5—  (ru  )  +  — 
dr  r  r 


9u„  9u 

-A  + _ z 

90  ^  9t 


(4-39) 


(4-40) 


The  o n e - d imen s i o n a  1  ax  i-symmet r ic  or  "pipe-flow" 
approximation  is  obtained  by  satisfying  the  continuity  and 
momentum  equations  on  an  average  by  integrating  them  through  the 
thickness.  The  resulting  one-dimensional  forms  of  the  continuity 
and  momentum  equations  are 


c2 


P 

9t  P 


2[b  u^(h)  -  a  uj.(a)] 


9u: 


(b^  - 


a^) 


=  0 


9x 


(4-41) 
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and 


,  5 

^  ^  _ X 

p  3x  "  3t 
rb 

prdr 


where 


2  2 
(b  -a^) 


(4-42) 


(4-43) 


^^rdr 


/t-2  2. 

(b  -  a  ) 


Upon  taking  the  temporal  derivative  of  Eq.  (4-41)  and  the  spatial 
derivative  of  Eq  (4-42),  and  combining  them,  we  are  led  to  the 
"pipe-flow"  equation 


(4-44) 


J_  „  2[b|^(b)  -a^(a)l 

3t2"9x2  > 


(4-45) 


For  the  axisymmetric  problem  presented  in  the  previous  section. 


(n=o) , 


U  (a)  =  io)  cos  ^ 


Uj.(b)  =  0 


(4-46) 


and  the  velocities  at  the  end  closures  vanish. 


u*(o)  =  u*(£)  =  0 

X  X 


(4-47) 


Thus,  the  boundary-value  problem  reduces  to 


^2  *  2 


p  *  .  u  -*  2aaj  p 
P  =  — ^ - 


■^^2  '  ^.2 
dX  C 


(b  -  a  ) 


,  mTTX  iut 

C04  e 


and 


3x  ■  |5-  “ 


(4-48) 
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and  its  solution 

is 

2 

Zao)  p  nrrrx  icjt 

P  ■  2  ® 

(b^  -  a^)[^  - 
c 

(4-49) 
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V.  Numerical  Calculations  and  Results 

The  geometry  of  the  capped  concentric  cylindrical  shell  used  in 
the  calculationa  is  shown  in  Figure  3.  For  illustrative  purposes 
three  types  of  internal  harmonic  excitations  have  been  considered: 

(1)  an  ax i s ymme t r i c  unit  pressure  acting  on  the  central  region  over 
a  length  of  22.5”,  (2)  a  unit  line  load  acting  along  the  central 

t 

region,  and  (3)  concentrated  axial  loads  acting  at  the  poles  of 
the  spherical  closures. 

The  fluid  pressure  field  has  been  obtained  for  nondimens ional 
frequencies  ranging  from  ka  =  0.18  to  3.52  for  the  single  shell  and  for  the 

shell  with  and  without  inner  fluid. 

Because  of  the  longitudinal  symmetry  of  the  loading  functions, 
only  longitudinally  symmetrical  structural  modes  and  surface 
expansion  functions  were  used.  Calculations  were  carried  out  by  using 
a  total  of  80  structural  modes  and  35  surface  expansion  functions  for 
each  of  the  first  four  circumferential  modes  (i.e.  for  n=0,1,2,3). 

The  impedance  functions  for  each  of  the  surface  expansion  functions 
were  obtained  by  subdividing  each  of  the  wetted  surfaces  into  180 
discrete  bands  (15  bands  for  each  sphere,  and  the  remaining  150 
bands  located  on  the  cylinder),  so  that  the  band  widths  were 
approximately  1.5  in. 

The  pressure  and  velocity  influence  coefficients  were  calculated 
by  using  surface  integrals  for  the  self-influence  coefficients,  and 
line  integrals  elsewhere.  However  for  the  axisymmetric  mode,  n=0,  all 
Influence  coefficients  pertaining  to  the  inner  fluid  were  evaluated  by 
means  of  surface  Integrals.  This  was  necessary  because  of  the 
extreme  sensitivity  of  the  axisymmetric  response  of  the  inner  fluid 
to  numerical  inaccuracies,  especially  at  low  frequencies. 
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All  calculations  wera  performed  on  the  CYBER  176.  The  execution 
time  required  to  obtain  80  structural  modes  for  a  particular 
circumferential  mode  number  was  approximately  30  secs  for  the  single 
shell  case/  and  90  secs  for  the  double  shell.  For  the  single  shell 
about  3  minutes  of  execution  time  was  required  to  calculate  the 
influence  coefficients  and  impedance  functions  associated  with  all  of 
the  surface  expansion  functions,  and  an  additional  20  sec  were 
required  to  calculate  the  structural  responses  and  the  far-field 
pressures.  These’ were  increased  to  20  minutes  and  40  secs, 
respectively  for  the  case  of  the  double  shell. 

The  physical  constants  used  for  all  calculations  are: 

FLUID. 

(Density)  P  =  9 . 45216xl0"^lb-sec2/in  ^ 

(Acoustic  Velocity) c  =  60,000  in/sec 
STEEL  SHELL 

(Density)  P=  73.3  xlO  ^Ib-sec^/in^ 

(Acoustic  Velocity)  c=  202,284  in/sec 

Results  for  the  far  field  pressure  at  points  located  32  radius 
(R  =  32a£)  from  the  geometric  center  of  the  structure  are  shown  in 

Figures  (5-21) 

The  radiated  pressure  field  for  the  axisymmetric  loading  of  the 
central  compartment  is  shown  in  Figs  (5-10)  .  The  results  indicate  that 
at  the  lower  frequencies  the  coupling  affected  by  the  inner  fluid  has 
little  effect  upon  the  radiated  pressure.  However,  no  general  trend 
can  be  observed  as  to  the  efficacy  of  using  a  particular  shell 
construction  to  reduce  the  far  field  radiation.  Indeed,  one  can  note 
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that  the  relative  radiation  characteristics  of  the  shell 
configurations  vary  with  frequency.  Of  course,  it  is  well  known  that 
the  relationship  between  the  forcing  frequency  and  the  in-fluid 
resonant  frequency  of  the  structure  is  the  primary  factor  in 
determining  the  magnitude  of  the  far  field  pressures. 

Figs  (11-15)  show  the  results  for  the  line  load  acting  on  the 
central  compartment,  while  Figs.  (16-21)  present  the  results  for  the 
end  cap  loadings.  For  these  cases,  as  for  the  previous  case,  no 
general  trend  is  observed. 

Because  of  numerical  difficulties  encountered  when  applying 
the  integral  methods  to  problems  in  which  internal  fluids  are 
present,  an  investigation  of  some  simple  approximations  to  the 
internal  fluid  problem  was  undertaken.  As  a  vehicle  for  this 
preliminary  investigation,  the  problem  of  two  concentric  circular 
cylindrical  shells  separated  by  a  small  gap,  which  is  filled  with 
an  acoustic  fluid,  was  considered  (Fig.  4).  The  outer  cylinder 
and  the  end  closure  rings  were  assumed  to  be  rigid,  and  the 
normal  displacement  of  the  inner  cylinder  was  assumed  to  vary 
cosinusoidally  in  both  the  circumferential  direction  and  along 
its  generator. 

Two  approximations  were  considered  —  a  first  order 
approximation,  in  which  the  fluid  potential  is  expressed  in  terms  of  a 
three-term  series  expansion  through  the  thickness,  and  a  simple  one¬ 
dimensional  or  "pipe  flow"  approximation.  A  comparison  of  the  results 
obtained  from  these  approximations  to  those  determined  analytically 
from  the  exact  formulation  is  given  in  Fig  (22  -  32).  It  should 
be  noted  that  the  very  large  pressures  obtained  are  simply  due 
to  the  fact  that  a  unit  displacement  amplitude  of  the  inner 
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cylinder  was  used. 

For  the  axisymmetric  case  (n=0),  Figs  (22-26),  the  first-order 
approximation  yields  excellent  results  over  a  wide  frequency  band  . 

The  lowest  frequency  for  which  the  approximation  is 

valid  varies  from  ka  =  0.088  at  the  axial  wave  number,  m  =  0  to  ka  =  1.58  at 
m=8  .  At  the  highest  frequency  f.or  which  calculations  have  been 
performed  (ka  =  4.40)  the  maximum  error  is  about  3%.  It  should  be  noted 
that  the  one-dimensional  approximation  yields  good  results  over  a  more 
restricted  frequency  range,  i.e.,  the  error  is  within  2%  for  ka  =  0.88 
to  l./b  tor  m=0  and  for  ka  =  1.67  to  2.11  for  m=8. 

For  the  asymmetric  cases  the  one-dimensional  approximations  is 
invalid.  However  the  first-order  approximation  yields  results  that 
are  within  2%  in  error  over  the  range  of  ka  =  0.088  -  4.40  for  n=l,  m=0  and 
over  the  range  of  ka  =  1.67  -  4.40  for  n=l,  m=8.  Similar  results  are 
obtained  for  the  n=2  modes. 

It  should  be  noted  that  the  first-order  approximation  yields 
excellent  results  over  the  frquency  range  where  the  exact 
solultion  is  expressed  in  terms  of  the  ordinary  Bessel  function, 
i.e.  when  toW  2  •  When  these 

functions  are  replaced  by  the  modified  Bessel  functions,  and  the 
three-term  series  approximation  no  longer  yields  meaningful 
results,  pointing  out  that  a  higher-order  series  approximation  is 
required . 

Thus  this  preliminary  investigation  of  the  first-order 
approximations  has  illustrated  that  such  a  series  expansion 
through  the  thickness  does  lead  to  a  significant  reduction  of  the 
effort  required  to  determine  the  far  field  radiation  characteristics 


-30- 


for  structures  which  contain  an  internal  fluid.  To  accomodate  a 
larger  frequency  rangOi  and  consider  more  complex  geometries  for 
which  an  internal  fluid  is  present,  a  higher-order  series 
expansion  is  required. 
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FIG.  16  RADIATED  FLUID  PRESSURE  (dB)  AT  R  *  32 

CONCENTRATED  LOAD  ON  ENCLOSURES  ;  f  »  100  Hz 


PRESSURE  IN  DECIBELS 


A«0* 

40  50  60  70  80  90 

PRESSURE  IN  DECIBELS 


45* 


50* 


□  SINGLE  SHELL 

O  DOUBLE  SHELL  (NO  INNER  FLUID) 
•  DOUBLE  SHELL  (INNER  FLUID) 


1 

00 

I 


FIG.  17  RADIATED  FLUID  PRESSURE  (dB)  AT  R  «  32 
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FIG.  18  RADIATED  FLUID  PRESSURE  (dB)  AT  R  »  32 
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FIG.  21  RADIATED  FLUID  PRESSURE  (dB)  AT  R  «  32 
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